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In this work we study interacting spinless fermions on a two-chain ladder with inter-chain pair 
tunneling while single-particle tunneling is suppressed at low energy. The model embodies a Z2 
symmetry associated with the fermion parity on each chain. We find that when the system is driven 
to the strong-coupling phase by the pair tunneling, Majorana excitations appear on the boundary. 
Such Majorana edge states correspond to two-fold degeneracy of ground states distinguished by 
different fermion parity on each chain, thus representing a generalization of one-dimensional topo- 
logical superconductors. We also characterize the stability of the ground state degeneracy against 
local perturbations. Lattice fermion models realizing such effective field theory are discussed. 

PACS numbers: 71.10.Pm, 03.67.Lx 



I. INTRODUCTION 

One-dimensional Topological Superconductors (TSC) 
are novel quantum phases of matter characterized by 
zero-energy Majorana edge states^. A lot of interest on 
TSC has been aroused due to the prospect of observ- 
ing Majorana particles in condensed matter systems^ as 
well as exploiting them as the building blocks of topo- 
logical quantum computers^! Various proposals of real- 
izing non-Abelian TSC in solid state systems have been 
put forward, e.g., in semiconductor/superconductor het- 
ero struc ture^HSl, anc l non-centrosymmetric superconduc- 
ton!™. 

The theoretical description of Majorana fermions in 
TSC is usually based on BCS mean-field Hamilto- 
nian which is essentially a non-interacting theory. An- 
other common feature of most existing studies of one- 
dimensional TSC is that the BCS pairing (thus the long- 
range superconducting order) is introduced by proxim- 
ity effect, since the strong quantum fluctuation in one 
dimension prevents spontaneous breaking of any contin- 
uous symmetry (the Mermin- Wagner theorem). The in- 
terplay between Majorana physics and interaction effects 
remains largely unexplored. Important questions such as 
how interactions affect Majorana fermions in TSC, and 
whether one-dimensional TSC can be induced from short- 
range interactions, have not been fully addressed. Sev- 
eral theoretical studies on the effects of interactions on 
Majorana fermions in proximity-induced TSC have been 
performed recently^"^, confirming the stability of Ma- 
jorana fermions against weak and moderate interactions. 
On the other hand, it is quite remarkable that the topo- 
logical classification of one-dimensional non-interacting 
fermionic systems with time-reversal symmetry is dra- 
matically changed by interactions^^. 

In this work we present a generic field-theoretical 
model of spinless fermions on two-chain ladders moti- 
vated by the second question that whether short-range in- 
teractions can induce TSC in one dimension. The model 



generalizes the simplest one-dimensional TSC, namely 
spinless fermions with p-wave pairing (also known as Ma- 
jorana chain)MJ, to interacting two-chain systems. Instead 
of introducing pairing by proximity effect, the effective 
field theory includes inter-chain pair tunneling with inter- 
chain single-particle tunneling being suppressed. There- 
fore the fermion parity on each chain is conserved. When 
the pair-tunneling interaction drives the system to strong 
coupling, localized Majorana zero-energy states are found 
on the boundaries, which represents a nontrivial many- 
body collective state of the underlying fermions. We then 
demonstrate that in a finite-size system the Majorana 
edge states lead to (nearly) degenerate ground states with 
different fermion parity on each chain, thus revealing its 
analogy with the Majorana edge states in non-interacting 
TSC. The degeneracy is shown to be robust to any weak 
intra-chain perturbations, but inter-chain single-particle 
tunneling and backscattering can possibly lift the degen- 
eracy. We also discuss a lattice model where such field 
theory is realized at low energy. 

II. FIELD-THEORETICAL MODEL 

We start from an effective field-theoretical description 
of the model for the purpose of elucidating the nature of 
the Majorana edge states. We label the two chains by 
a = 1,2. The low-energy sector of spinless fermions on 
each chain is well captured by two chiral Dirac fermions 
^ 'l I R.a( x ) ■ The non-interacting part of the Hamiltonian 
is simply given by Hq = J dx Ho (x) where 

Ho = -Wp ^ [^Ra^Ra ~ ^L a d xi>La) ■ (1) 
a 

Four-fermion interactions can be categorized as intra- 
chain and inter-chain interactions. Intra-chain scattering 
processes (e.g., forward and backward scattering) are in- 
corporated into the Luttinger liquid description of spin- 
less fermions and their effects on the low-energy physics 
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are completely parameterized by the renormalized veloci- 
ties v a and the Luttinger parameters K a . We assume that 
the filling of the system is incommensurate so Umklapp 
scattering can be neglected. For simplicity we assume the 
two chains are identical so v\ — V2 — v, K\ — K2 = K. 

We now turn to inter-chain interactions. Those that 
can be expressed in terms of the densities of the chiral 
fermions can be absorbed into the Gaussian part of the 
bosonic theory after a proper change of variables(see be- 
low) and we do not get into the details here. We have to 
consider the pair tunneling and the inter-chain backseat - 
tering: 



Tibs = gh S {i> L1 i>Ri4> R 2i>L2 + 1 «-> 2). 



(2) 



The microscopic origin of such terms is highly model- 
dependent which will be discussed later. The motivation 
of studying pair tunneling is to "mimic" the BCS pairing 
of spinlcss fermions without explicitly introducing super- 
conducting pairing order parameter. 

The Hamiltonian of the effective theory is then ex- 
pressed as H = Hq + Whs + "Hpair- Notice that total 
fermion number N = Ni + N2 is conserved by the Hamil- 
tonian, but Ni and N2 themselves fluctuate due to the 
tunneling of pairs. However, their parities (— l)^ a are 
still separately conserved. Due to the constraint that 
(-1)* 1 ■ (-1)* 2 = (-l) N , we are left with an overall 
Z2 symmetry. Therefore we define the fermion parities 
P a = (— 1) °, the conservation of which is crucial for 
establishing the existence and stability of the Majorana 
edge states and ground state degeneracy. In the following 
we refer to this overall Z2 fermion parity as single-chain 
fermion parity. It is important to notice that the con- 
servation of the single-chain fermion parity relies on the 
fact that there is no inter-chain single-particle tunneling 
in our Hamiltonian. We will address how this is possible 
when turning to the discussion of lattice models. 

We use bosonizatior l 18 * 19 ! to study the low-energy 
physics of the model. The standard Abelian bosoniza- 
tion reads 



\/27rao 



(3) 



where ao is the short-distance cutoff, r = +/— for R/L 
movers and r) r ^ a are Majorana operators which keep track 
of the anti-commuting character of the fermionic opera- 
tors. We follow the constructive bosonization as being 
thoroughly reviewed in [T9]. The two bosonic fields f a 
and 9 a satisfy the canonical commutation relation: 



[d x ipa(x),0 a (x')] =lS(x-x'). 



(4) 



The if a field is related to the charge density on chain a 
by Pa — -^d x ip a , and 9 a is its conjugate field, which can 
be interpreted as the phase of the pair field. 



It is convenient to work in the bonding and anti- 
bonding basis: 

<f± = ^(<Pi±iP2),e ± = -^(6i±9 2 ). (5) 

The resulting bosonized Hamiltonian decouples as % — 
U++U-: 



H- 



V -±[K + (d x 6 + f + K+ 1 (d x v + ) 2 }, 

V — [K-idJ-Y + KZ 1 ^-?] 

—, gp .„ cos %/8~7r#_ + - , gbs .„ cos V8nip. 



(6) 



Here ao is the short-distance cutoff. This decoupling of 
the bonding and the anti-bonding degrees of freedom is 
analogous to the spin-charge separation of electrons in 
one dimension. Without any inter-chain forward scatter- 
ing, we have K± — K, v± = v. 

The bonding sector is simply a theory of free bosons. 
The Hamiltonian in the anti-bonding sector can be ana- 
lyzed by the perturbative Renormalization Group (RG) 
method, assuming the bare couplings g p and gbs are 
weak. RG flow of the coupling constants are governed 
by the standard Kosterlitz-Thouless equations^ 



dyp 
dl 

dj/bs _ 
d/ 
din if. 

dl 



(2 - 2KI 1 )y p 
-- (2-2X_)t/ bs 



(7) 



2KZV-, 



where w_ = , Vh s — are the dimensionless cou- 

& 7tv— ' & u& Tri- 
pling constants and I — In ^ is the flow parameter. 

When K- > 1 (corresponding to attractive intra-chain 
interaction), y p is relevant and flows to strong-coupling 
under RG flow, indicating gap formation in the anti- 
bonding sector, while j/bs is irrelevant so can be neglected 
when considering long-wavelength, low-energy physics. 
Semiclassically, the 9- is pinned in the ground state. 
From now on, we will assume K— > 1 and neglect the 
irrelevant coupling ?/bs- 



III. MAJORANA ZERO-ENERGY EDGE 
STATES 

To clarify the nature of the gapped phase in the anti- 
bonding sector, we study the model at a special point 
K- = 2, known as the Luther-Emery pointPH, where the 
sine-Gordon model is equivalent to free massive Dirac 
fermions. First we rescale the bosonic fields: 



<P- 



V- 



(8) 



and define the chiral fields by (p r ^ = \{(p- + r6~). 
Neglecting the irrelevant backscattering term, "H_ is 
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refermionized to 



iv-(x f R d x x R - x\9 x xl) + imiXRxl ~ XlXr) 



t ,->t 



where the Dirac fermionic fields Xr are given by 
1 



Xr 



y/2ira 



iry 47T(p r _ 



(9) 

(10) 



with the fermion mass m — £ r are again Majorana 
operators. It is quite clear that effective theory Q also 
describes the continuum limit of a Majorana chain, which 
is known to support Majorana edge states-. 

However, caution has to be taken here when dealing 
with open boundary condition(OBC). We impose open 
boundary condition at the level of underlying lattice 
fermionic operators^ 



' 1pLa(x) 



(ii) 



where c, a are annihilation operators of fermions and 
x = idQ. Since the chain terminates at x = and x = L, 
we demand £q — c^v+i = where N = L/ao is the num- 
ber of sites on each chain. Let us focus on the boundary 
x = 0. Thus the chiral fermionic fields have to satisfy 
"0i?.a(O) = — V'ia(O). Using the bosonization identity, we 
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ary condition of the anti-bonding field: 



find V'a(O) = from which we can deduce the bound- 



P-(0)=0. 



(12) 



Therefore, we obtain the boundary condition of the 
Luther- Emery fermionic fields as xb,(0) — Xl(0). The 
Hamiltonian is quadratic in \ and can be exactly diago- 
nalized by Bogoliubov transformation. We find that the 
Luther-Emery fields have the following representation: 



Xl(x) 



(13) 



Here . . . denotes the gapped quasiparticles whose forms 
are not of any interest to us. The 7 is a Majorana 
field(i.e., 7 = 7T) and because ["H_,7] = 0, it represents 
a zero-energy excitation on the boundary. 

Now suppose the system has finite size L> £ = V-/m. 
The same analysis implies that we would find two Ma- 
jorana fermions localized at x — and x — L respec- 
tively, denoted by 71 and 72. As in the case of TSC, the 
two Majorana modes have to be combined into a (nearly) 
zero-energy Dirac fermionic mode: c= ^(71+272). Oc- 
cupation of this mode gives rise to two degenerate ground 
states. Tunneling of quasiparticles causes a non-zero 
splitting of the ground state degeneracy: AE me _i ^ 

|23|24I 

We notice that very similar techni que was previously 
applied to the spin- 1/2 edge excitation d 22 | 25 | 2e | in the Hal- 
dane phase of spin-1 Heisenberg chain, the §Q(n) spinor 
edge states in the SO(n) spin chairP^ and also th e edg e 
state in an attractive one-dimensional electron ga d 28 * 29 l 



To understand the nature of the Majorana edge state, 
we have to explicitly keep track of the Klein factors which 
connect states with different fermion numbers. Therefore 
we separate out the so-called zero mode in the bosonic 



field 



7 X T: 



where 



and write ij) ra 

the Klein factors F ra are bosonic operators that decrease 
the numbers of r-moving fermions on chain a by one^l. 
Bosonized form of ^ has a product of the Klein factors 
FIjoF^FijiFri in it. Since this term is to be refermion- 
ized as ~ XlXr-i we are naturally led to define new Klein 
factors F r = Fj 2 F r i for \r- Notice that so-defined Klein 
factors satisfy {P a ,F r } = 0, i.e. F r change single-chain 
fermion parity. Then the fermionic fields that refermion- 
ize the sine-Gordon theory at the Luther-Emery point 
should take the form 



1 



£ p irs/4nip r 



(14) 



Thus one can identify that Xr corresponds to inter-chain 
single-particle tunneling. The ground state \G) of the 
Hamiltonian ^ can be schematically expressed as 



\G) = exp 



dxxdx2 X* (xi)g(xi, x 2 )xH x ^) 



vac) 



(15) 

where g{x\,X2) is the Cooper-pair wave function of the 
spinless p-wave superconductor and |vac) is the vacuum 
state of x fermion. With the definition (14), it is easy 



to check that \G) is a coherent superposition of Fock 
states having the same single-chain fermion parity, thus 
an eigenstate of P a . On the other hand, the Majorana 
fermion 7, being a superposition of x aud X* > changes the 
single-chain fermion parity: {7, P a } = 0. As a result, the 
two degenerate ground states \G) and $\G) have different 
single-chain fermion parity which is the essence of the 
Majorana edge states. If the total number of fermions N 
is even, then the two (nearly) degenerate ground states 
correspond to even and odd number of fermions on each 
chain, respectively. 

So far all the conclusions are drawn at the Luther- 
Emery point K_ = 2. Once we move away from the 
Luther-Emery point, the theory is no longer equivalent 
to free massive fermions. An intuitive way to think about 
the situation is that if we move away from the Luther- 
Emery point, the x fermions start to interact with each 
other. Since the Majorana edge states are protected by 
the bulk gap as well as the single-chain fermion par- 
ityEHH] we expect the qualitative features hold for the 
whole regime K_ > 1 based on adiabatic continuity. 

Notice that the bonding sector remains gapless. In our 
field-theoretical model, the bonding and anti-bonding de- 
grees of freedom are completely decoupled so the gapless- 
ness of the bonding boson does not affect the degeneracy 
in the anti-bonding sector. 
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FIG. 1: Schematic view of the two chains coupled by pair 
tunneling(denoted by dashed lines). The chains are bended 
near the two ends to avoid the single-particle tunneling. 



IV. STABILITY OF THE DEGENERACY 

We now examine whether the ground state degener- 
acy we have found has a topological nature. Here we 
define a topological degeneracy of the ground states by 
the following criteria: the two degenerate ground states 
are not distinguishable by any local order parameters(i.e. 
the difference of the expectation values of any local order 
parameters in the two ground state must be exponen- 
tially small in system size ). By local, we mean local 
operators in the original fermionic operators , otherwise 
we can easily find such an operator in the bosonic rep- 
resentation. For example, in the model ^ the operator 
O(x) — cos y/2n9-(x) can distinguish the two degenerate 
ground states. But the operator itself is highly non-local 
in terms of the original fermionic operators. 

First of all, by analogy with Majorana chain it is quite 
obvious that any local operators that involve even num- 
bers of fermion operators on each chain are not able to 
distinguish the two ground states because such operators 
always commute with single-chain fermion parity opera- 
tor. Therefore we only have to consider operators that 
consist of odd number of single-chain fermion operators. 
They change the single-chain fermion parity and thus 
presumably connect the two degenerate ground states. 
Since all such operators can be decomposed into prod- 
ucts of single-particle inter-chain tunneling and backseat - 
tering operators, it is sufficient to consider these single- 
particle operators. 

Let us start with single-particle inter-chain tunneling 

Ot= £ ($ P & r + h.c). (16) 

r=R,L 

Its bosonic representation is 

2 

£>t = cos \/2^_ cos V2tt6_ . (17) 

7ra 

First let us consider the case when the operator is taken 
in the bulk of the chain away from any of the boundaries. 



Because 9- is pinned in the ground states, tp- gets to- 
tally disordered and therefore (Ot) oc (cos V^tt^P-) = 0, 
which is just equivalent to the fact that the Luther-Emery 
fermions are gapped. However, this is no longer true as 
one approaches the ends of the chains, since there exists 
zero-energy edge states. Let us focus on the left bound- 
ary x = 0. The boundary condition of the anti-bonding 
boson field ip _ has been derived: <fi-(0) = 0. With the 
boundary condition, we proceed with Luther-Emery so- 
lution at K- = 2 and find O T (0) ~ x(0) + x f (0). Thus 
Ot (0) has nonvanishing matrix element between the two 
ground states, independent of the system size. As a re- 
sult, the two-fold degeneracy is splitted. 

We now turn to the inter-chain backscattering 

Cb = i'ln^PiL + i'l^iR + h- c - 

2 >_ >_ (18) 

= cos v 2n(p + cos v 27r6L . 

7ra 

An analysis similar to the single-particle tunneling leads 
to the conclusion that backscattering at the ends also 
splits the degeneracy. However, even if the backscatter- 
ing occurs in the middle of the chain, it still causes a 
splitting of the ground states decaying as a power law 
in system size L. To see this, let us consider a sin- 
gle impurity near the middle of the chain, modeled by 
Ob(x) where x s» L/2. We assume that the backscat- 
tering potential is irrelevant under RG flow and study 
its consequence. The splitting is then proportional to 
(cos y/2~TTip + (x)) since cos \/2tt9- has different expecta- 
tion values on the two ground states. Because ip+ is 
pinned at x = 0, (cos \/2ttip + (x)) ~ l/x K +. Therefore 
the splitting of the ground states due to a single impu- 
rity in the middle of the system scales as 1/ L K + . 

We thereby conclude that the ground state degeneracy 
is spoiled by the single-particle intcr-chain tunneling near 
the boundaries and the backscattering processes in the 
bulk. To avoid the unwanted tunneling processes near 
the ends, one can put strong tunneling barriers between 
the two chains near the ends, or the chains can be bended 
outwards so that the two ends are kept far apartP^, as 
depicted in Fig. [I] 

V. LATTICE MODEL 

We now show that the field theory ^ can be realized 
in lattice models of fermions. We consider th e mo del 
of two weakly coupled chains of spinless fermionj 18 * 31 * ^. 
The Hamiltonian reads 

H = -t^(cl +l a c ia +h.c.) + £ V(r)h la n i+r . a 

i,a i,a,r ,^ n x 

t ^ ' 

-t x 2^(q 2 Cii + h-c)- 

i 

Here a = 1, 2 labels the two chains. We assume the filling 
is incommensurate to avoid complications from Umklapp 
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scatterings. V(r) is an intra-chain short-range attrac- 
tive interaction between two fermions at a distance r (in 
units of lattice spacing). Thus without inter-chain cou- 
pling, each chain admits a Luttinger liquid description 
with two control parameters: charge velocity v and Lut- 
tinger parameter K (we assume V is not strong enough 
to drive the chain to phase separation). 

We bosonize the full Hamiltonian and write the theory 
in the bonding and anti-bonding basis. Hamiltonian in 
the bonding sector is just a theory of free bosons. In the 
anti-bonding sector, it reads 



H 



K(d x 9) 2 + ^-(d x ^) 



2*j 
na 



COS V 27T(/J cos v 27:9. 

(20) 

The bosonic fields (p and are in the anti-bonding basis. 
The perturbation (t±) term has nonzero conformal spin 
which implies that two-particle processes are automati- 
cally generated by RG flow even when they are absent 
in the bare Hamiltonian. Therefore, one has to include 
two-particle perturbations in the RG flow 



n 2 



gi 



(7ra ) : 



■ cos v 8ir<p 



fj2 



(7ra ) : 



cos v8tt8 



(21) 



The RG flow equations for weak couplings have been de- 
rived by Yakovenk o |3"T] and Nersesyan et al. [35]. Here 
we cite their results^ 



dz 
dl 
dj/i 

d; 

dj/2 

d; 

dK 
Id 



K + K- 



(2 - 2K)y 1 + {K - K~ l )z 2 
(2-2K- 1 )y 2 + (K- 1 -K)i 
\ivl-vlK 2 ) 



(22) 



where the dimensionless couplings are defined as z — 
and y 12 = Shl. 

Since we are interested in the phase where the pair 
tunneling dominates at low energy, we assume K > 1 
so yi is irrelevant and can be put to 0. Also we neglect 
renormalization of K. Integrating the RG flow equations 
with initial conditions z(0) = Zq ^ 1, y 2 (0) =0we obtain 



2/2(0 



,K~ 



K 



2a 



e 2{l-a)l_ e 2{l-K-i)U (23) 



where a = h(K + K 1 — 2). Assume K 1 < a, then 

the large-? behavior of y 2 is dominated by e 2< - 1 ~ K y 2 
becomes of order of 1 at I* sa — lnzo/(l — K^ 1 ), where 

the flow of z yields z(l*) « z^ a K K ' < 1 given 

ZD < 1. This means that if K > y/2 + 1 (so K~ l < a), 
then y 2 reaches strong-coupling first. Thus the strong- 
coupling field theory is given by (|6|. 



VI. CONCLUSIONS 



To conclude, we consider the strong-coupling phases of 
model of spinless fermions on a two-chain ladder driven 
by the pair tunneling. We find that, through Luther- 
Emery solution of the strong-coupling model, there ex- 
ists zero-energy excitations on the edges of the ladder 
represented by Majorana fermions. On a finite system 
there are always two such Majorana edge states which 
can be combined to a Dirac fermionic mode and therefore 
the ground states are two-fold degenerate, correspond- 
ing to the mode being occupied or unoccupied. We fur- 
ther clarify the nature of the ground state degeneracy 
and show that the two states have different fermion par- 
ity on each chain. This is in complete analogy with the 
one-dimensional topological superconductor, where there 
are two ground states with different total fermion parity. 
However, in our case, the one-dimensionality prevents the 
spontaneous breaking of the global U(l) symmetry and 
what we find is the degeneracy between the states with 
different fermion parity on each chain, subject to the con- 
straint that the total number of fermions is fixed. This 
is an important distinction between the strong-coupling 
phase studied in this work and the one-dimensional topo- 
logical superconductor. What is more, the degeneracy we 
have found is a purely interaction effect and thus goes be- 
yond the mean-field theory of topological superconduc- 
tivity in one dimension (essentially non-interacting). 

We further characterize the robustness of the ground 
state degeneracy. We find that the degeneracy is im- 
mune to any local perturbations that preserve the single- 
chain fermion parity. The inter-chain single-particle tun- 
neling in the bulk is prohibited by the existence of a 
single-particle gap as well. However, near the bound- 
aries the bulk gap vanishes (hence the existence of zero- 
energy states) and inter-chain single-particle tunneling or 
backscattering can lift the degeneracy by a finite amount 
that is independent of the system size. Furthermore, due 
to the gaplessness of the bonding sector, the inter-chain 
backscattering in the bulk also change the splitting of the 
degeneracy to be power-law in system size. 

We also discuss a lattice model of two weakly coupled 
spinless fermion chain where such low-energy effective 
field theory is realized. We show that there is a range 
of the Luttinger parameter K such that the inter-chain 
single-particle tunneling becomes irrelevant (or less rel- 
evant) under RG flow, but the two-particle pair tunnel- 
ing, generated by the single-particle tunneling, becomes 
relevant and grows to strong coupling eventually. This 
confirms the validity of our general field-theoretical ap- 
proach. 

Note added. During the finalization of the manuscript, 
we learnt that related works in the context of spin-orbit 
coupled nanowires has been done by Fidkowski et al^^ 
and also by Sau et. alW^ . 
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